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Radio- quiet quasar environments at 0.5 < z < 0.8 
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ABSTRACT 

We have quantified the galaxy environments around a sample of 0.5 < z < 0.8 radio- 
quiet quasars using the amplitude of the spatial galaxy-quasar correlation function, 
Bgq. The quasars exist in a wide variety of environments, some sources are located 
in clusters as rich as Abell class 1-2 clusters, whereas others exist in environments 
comparable to the field. We find that on average, the quasars prefer poorer clusters 
of w Abell class 0, which suggests that quasars are biased tracers of mass compared 
to galaxies. The mean i3gq for the sample is found to be indistinguishable from the 
mean amplitude for a sample of radio-loud quasars matched in redshift and optical 
luminosity. These observations are consistent with recent studies of the hosts of radio- 
quiet quasars at low to intermediate redshifts, and suggest that the mechanism for the 
production of powerful radio jets in radio-loud quasars is controlled by processes deep 
within the active galactic nucleus itself, and is unrelated to the nature of the hosts or 
their environments. 
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1 INTRODUCTION 

The study and interpretation of the differences and simi- 
larities between radio-loud and radio-quiet quasars (RLQs 
and RQQs, hereafter) has kept astronomers busy for a long 
time. The spectral energy distributions of RLQs and RQQs 
differ only in points of detail at all wavelength bands, apart 
from the radio. Both quasar types possess broad emission 
lines and have a stellar appearance in even the highest res- 
olution optical images. They also show the same evolution 
in number density with a peak at 2 ~ 2-3, but the RQQs 
outnumber the RLQs by a factor of 10. RQQs have radio 
luminosities 2-3 orders of magnitude lower than their radio- 
loud counterparts, and on radio maps they are compact with 
typically only a weak, flat-spectrum component coincident 
with the optical nucleus. On the other hand, the RLQs have 
extended lobes of radio emission with hotspots at the outer 
edges of the radio structure. The radio-emitting lobes are 
fed by powerful jets emerging from a bright, central core, 
but RQQs can also have jet-like structures (Blundell et al. 
1996; BlundeU & Beasley 1998), ahhough with bulk kinetic 
powers ~ 10'^ times lower than for RLQs (Miller, Rawlings & 



SEmail: wold@astro.su.se 



Saunders 1993). This suggests that both quasar populations 
have jet-producing central engines, but that the efficiency 
of the jet production mechanism is very different in the two 
cases. 

One way to learn more about how the two quasar types 
relate, is to study their host galaxies and their galactic en- 
vironment. A longstanding belief that RQQs are hosted by 
spiral galaxies and RLQs by ellipticals is now being ques- 
tioned. Recent studies (e.g. Dunlop et al. 1993; McLeod & 
Rieke 1994; Taylor et al. 1996; BahcaU et al. 1997; Boyce et 
al. 1998; McLure et al. 1999; Hughes et al. 2000) have found 
that powerful quasars at z ^ 0.5, both radio-loud and radio- 
quiet, seem to exist in galaxies with luminosities > L* , the 
luminosity at the break in the luminosity function. Still, a 
clear picture has not emerged, with some studies claiming a 
high fraction of disk morphologies amongst the radio-quiets 
(e.g. Percival et al. 2000), whilst others suggest that nearly 
all quasars are in giant ellipticals (e.g. McLure et al. 2000). 
Besides the technical difficulties associated with studying 
host galaxies, quasar optical luminosity and selection ef- 
fects may affect the results. E.g. the perhaps somewhat ar- 
tificial distinction between Seyfert galaxies and quasars at 
Mv = —23 leads to confusion as to whether an active galac- 
tic nucleus (AGN) is a quasar or a Seyfert. Seyfert galaxies 
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are spirals, are found at lower redshifts and have lower AGN 
luminosities than RLQs, so very few RLQs exist for compar- 
ison. Selection effects may therefore have caused spirals to 
be preferred as hosts of RQQs and giant ellipticals to be 
preferred as hosts of RLQs. 

Quasar environments have also been investigated on 
much larger scales than those of typical host galaxies. RLQs 
at 0.5 ^ z ^ 0.8 are often found to be associated with 
groups or clusters of galaxies (Yee & Green 1987; Ellingson, 
Yee & Green 1991a; Wold et al. 2000 (hereafter Paper I)), 
with richnesses varying from poor groups to Abell Class 1 
clusters or richer. This is perhaps not surprising since gi- 
ant elliptical galaxies, found to be the hosts of powerful 
RLQs, frequently reside in the centres of galaxy clusters. 
But what is the galaxy environment like for RQQs in this 
redshift range, and how does it compare with the environ- 
ment around RLQs with comparable AGN luminosities? 

In this paper, we try to answer this question by inves- 
tigating a sample of RQQs matched in redshift and optical 
luminosity to a sample of RLQs analysed in Paper L We 
quantify the galaxy environment using the amplitude of the 
spatial galaxy-quasar cross-correlation function, Bgq, often 
referred to as the 'clustering amplitude'. The first attempt 
at comparing the environments around RLQs and RQQs 
using the clustering amplitude was made by Yee & Green 
(1984), who found only a marginal difference. Later, Yee & 



Green (1987) did deeper imaging of the same RQQ fields, 
and added more RLQ fields to the sample. They repeated the 
analysis with improved parameters, and obtained an even 
smaller difference, but due to the small number of RQQs 
(seven) in their sample, they were unable to draw any firm 
conclusions. 

Th e problem was later addressed by Ellingson et al. 



(1991a), who also used the clustering amplitude to quan- 
tify the Mpc-scale environment around a sample of 65 
17.5 < m < 19 quasars at 0.3 < z < 0.6, of which half were 
radio-quiet. They found that the RLQs exist more often in 
rich galaxy environments than the RQQs do, which was be- 
lieved to support the RQQ/spiral host-RLQ/elliptical host 
hypothesis. Since this study, there has been little work aimed 
at comparing the environments of the two quasar popula- 
tions at moderate redshifts, but at both lower and higher 
redshifts several studies indicate that contrary to the con- 
clusion reached by Ellingson et al., the galaxy environment 
around RLQs and RQQs is similar. But also here, a clear pic- 
ture has not yet emerged, since several other studies seem 
to support Ellingson et al's conclusion. 

Some investigators use the clustering amplitude, which 
makes comparison between studies similar to ours easier, 
whereas other studies have other aims, and therefore use 
different techniques and methods. These studies can basi- 
cally be divided into two categories. The first category in- 
cludes measurements of the angular galaxy-quasar cross- 
correlation function by cross-correlating quasar catalogues 
with galaxy counts. The aim is to establish whether opti- 
cally and X-ray selected quasars (i.e. mostly RQQs) can be 
used as unbiased tracers of galaxies. This is important to 
know, since t he large quasars su rveys currently in progress, 
like the 2dF (IShanks et al. 200c| ) and the Sloan Digital Sky 
Survey (Fan et al. 1999), will be used to study large-scale 
structure out to high redshifts, but require knowledge about 
how quasars are biased with respect to galaxies. The sec- 



ond category includes more or less qualitative searches for 
galaxies in quasar and radio galaxy fields at both high and 
low redshifts. E.g. data obtained for host galaxy studies, at 
typically z ^ 0.3, can be utilized for environment studies, 
but the high resolution imaging required to study hosts of- 
ten results in a small field of view and therefore a relatively 
small sampling radius. 

Investigations of the galaxy-quasar angular correlation 
function are often concerned with scales larger than typi- 
cal clu ster scales. At low redshifts. Smith, Boyle & Maddox 
selected 169 {z < 0.3) optically identified quasars 



( [1995| ) 



from the EMSS (Einstein Medium Sensitivity Survey), and 
cross-correlated their positions with Bj < 20.5 galaxies from 
the APM galaxy catalogues in a region of radius 0?5 centered 
on the quasar. They found that the amplitude of the angu- 
lar galaxy-quasar cross-correlation function was identical to 
that of the APM galaxy angular correlation function, imply- 
ing that the quasars inhabit regions of space similar to that 
of normal galaxies. (At z — 0.3, 0?5 corresponds to alm ost 10 
Mpc). In a follow-up study. Smith, Boyle & Maddox ( |200C| ) 
cross-correlated X-ray selected quasars at 0.3 < z < 0.7 
with V < 23 galaxies on CCD images. They found a 2a 
excess around the z < 0.5 quasars on scales < 1 arcmin, 
but no excess for the z > 0.5 quasars. As an overall result, 
they claimed that the galaxy-quasar correlation function is 
consistent with that for faint galaxies. The quasars in their 
sample span a wide range in luminosity, from Mv ~ — 27 
to —20 {Ho = 50 kms"'^ Mpc~^), but no difference between 
high and low luminosity sub-samples was observed. 

Groom & Shanks (1999) cross-correlated 150 quasars 
selected both in the optical and in X-rays with bj < 
23 galaxies on AAT (Anglo- Australian Telescope) photo- 
graphic plates. They found the angular cross-correlation to 
be marginally negative, indicative of a small anti-correlation, 
possibly due to gravitational lensing. Most of the quasars in 
their sample lie at redshifts 1.0-2.2, and they claimed to 
reach M* at z « 1 and M* — 1 at z ~ 1.5. Possibly, this is 
not deep enough, as even M* might be too shallow into the 
l uminosity function to de tect excess galaxies ai z = 1.0-2.2 
(lYee & L6pez-Cruz 19991). 



Searches for galaxies in quasar and radio galaxy fields 
in connection with host s tudies have been conducted by e.g. 
Smith & Heckman (199C) who generated catalogues of com- 
panion galaxies within 100/i~^ Mpc of 31 z < 0.3 quasars 
and 35 powerful radio galaxies. They found that RQQs only 
have half the number of nearby, bright companions as RLQs 
and radio galaxies at the 95 per cent confidence level. Some- 
what inconsistent with Smith & Heckman's result is the 



result obtained by Fisher et al. (1996) who found no sig- 



nificant difference in the environments around optically lu- 
minous RLQs (six) and RQQs (14) at z ^ 0.3. Supporting 
Fisher et al's result is the recent study by McLure & Dunlop 
( [2000| ), who quantified the environments around 44 power- 
ful, z « 0.2 AGN, of which 13 were RLQs and 21 were 
RQQs (their sample includes Fisher et al's quasars). They 
found the mean clustering amplitude for the two quasar sam- 
ples to be consistent with each other, and corresponding to 
Abell richness class clusters. Dunlop et al. (1993), study- 



ing quasar hosts in the near-infrared, also discuss companion 
galaxies. Their sample consisted of 14 RLQs and 14 RQQs 
at z < 0.35 matched in redshift and V magnitude. They 
found that the surface density of bright {K > 18) compan- 
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ion galaxies was a factor of ~ 2.5 higher than in the field, 
and the excess was similar for both RLQs and RQQs. 

The high redshift studies of quasar environments typi- 
cally concentrate at z ^ 1. Quasars (and radio galaxies) are 
used as pointers to high redshift galaxies in order to learn 
about galaxy formation and evolution and to search for high 
redshift (proto)clusters. Galaxies, and even clusters, are cer- 
tainly found at 2 ^ 1 using this method (e.g. Dickinson 1994; 
Dehorn et al. 1997; HaU & Green 1998; Tanaka et al. 2000, 
Chapman, McCarthy & Persson 2000; Best 2000). 

Hi ntzen, Romanishin & Valdes (1991) and Boyle & 
Couch ( |l993 ) investigated 0.9 < 2 < 1.5 RLQs and RQQs, 
respectively. Hintzen et al. found a significant galaxy ex- 
cess within a 15 arcsec radius (~ 160 kpc) centered on 



the quasars in their sample, whereas Boyle & Couch ( 199J ) 
found no evidence for excess. However, Boyle & Couch's 
RQQs are 20 times fainter in the optical than the RLQs 
studied by Hintzen et al., and Boyle & Couch discuss this 
as a possible cause for the discrepancy. 
In contrast to this, Hutchings, Crampton & Johnson 



(1995) found on average a. > 5a excess of galaxies around 
both RLQs and RQQs (14 in total) at z ~ 1.1. They saw 
no significant difference between RLQ and RQQ environ- 
ments, which was found to consist of compa ct groups of 
starbursting galaxies. Also, Hutchings (1995a), who looked 
at galaxy companions to seven 2 = 2.3 quasars (only one 
RLQ) in the _R-band found an overall > 5a excess in a ~ 1 
arc min re gion around the quasars. However, Hutchings et 
al. (1999), using near-infrared imaging of both RLQs and 
RQQs at 0.9 < z < 4.2, note that the RQQs occupy poorer 
environments than the RLQs, but do not draw any firm 
conclusions due to the sma ll an d still incomplete sample. 
Teplitz, McLean & Malkan (1999) also report that the num- 
ber counts of galaxies in RQQ fields at 0.95 < z < 1.5 in 
J- and Tt'-band is comparable to the field counts, but their 
data is not very deep (only reaching ~ M* at 2 « 1). 

In summary, although the results seem to differ, the 
consensus about hosts and environments of RQQs - that 
they are typically late-type galaxies in poor environments 
- is being increasingly questioned. Several recent observa- 
tions at high and low redshifts seem to indicate that RQQs 
live in richer than average galaxy environment, and powerful 
RQQs at 2 ^ 0.5 are found in luminous, possibly early-type, 
galaxies. 

Apart from Ellingson et al's (1991a) study at 0.3 < 
z < 0.6, there have been no attempts to investigate what 
environments RQQs prefer at intermediate redshifts. We 
have therefore obtained data on a sample of 20 RQQs at 
0.5 < z < 0.8, and quantified the galaxy richness in the 
quasar fields using a similar analysis as Ellingson et al. The 
sample is described in Section |2|, where we also discuss the 
radio loudness of some of the quasars. In Section ^, we 
describe the observations and our observing strategy. The 
control fields used in the analysis are also commented upon 
here, especially how the galaxy counts vary with galactic lat- 
itude. The analysis using the clustering amplitude is briefly 
overviewed in Section % before the results are presented in 
Section |^. In Section we examine the colours and the 
radial distribution of the galaxies in the quasar fields. Using 
the clustering amplitude allows us to compare our results 
with some of the above mentioned studies, and this is done 
in Section ^, and in Section ^ we are able to directly and 



consistently compare with the sample of RLQs that was pre- 
sented in Paper I. Finally, we discuss the results in Section ^ 
and draw the conclusions in Section ^ 

Our assumed cosmology has Ho = 50 kms^^Mpc"^, 
fio = 1 and A = 0. 



2 THE SAMPLE 

The RQQs were selected from three optical surveys with 
different fiux limits in order to cover a wide range in AGN 
luminosity (B-band) within the given redshift range. Eight 
of the quasars are from the faint {B ^ 21) Durham/AAT 
UVX survey of Boyle et al. (|l990|) (designated 'BFSP') and 



ten are from the intermediate luminosity (16 ^ bj ^ 18.9) 
Large B right Quasar Survey (LBQS) by Hewett, Foltz & 
Chaffee (1995). There are also two high- luminosity quasars 



in the sample, selected from the B ^ 16.6 Bright Quasar 
Survey (BQS) by Schmidt & Green ([l983|). The sample is 



listed in Table 1 along with information about the observa- 
tions. 

The BQS and the Durham/AAT UVX survey are based 
on the UV-excess selection technique which finds z < 2.2 
quasars with stellar-like appearance and excess emission 
in the UV. The LBQS is an objective prism survey, from 
which quasar candidates are essentially selected by their blue 
colour. The LBQS and the Durham/AAT UVX survey are 
complete, but the BQS s urvey has been found to be incom- 
plete by a factor of 3.4 (Goldschmidt et al. 1992), possibly 
due to photometry errors. Also, the BQS is known to have an 
anomalously high radio-loud fraction at bright optical mag- 
nitudes, AIb < —24, compared to other optically selected 
samples (LaFranca et al. 1994), and Hooper et al. (199C) 
find that the radio-loud fraction at faint optical luminosities 
is lower in the BQS than in the LBQS. We have checked 
that inclusion of the two BQS quasars does not affect our 
results. 

The RLQ sample that we compare with in Sec- 
tion 1^ consists of 21 steep-spectrum RLQs at 0.5 < 
z < 0.8 covering the wide radio luminosity range 23.8 < 
log (L408MHz/WHz~^ sr^^) < 26.7. This sample was drawn 
from two different radio/optical flux-limited samples, the 
Molonglo/APM Quasar Survey (Serjeant 1996; Maddox et 
al., in preparation; Serjeant et al., i n preparation) and the 
7C Quasar survey (Riley et al. 1999). For details about the 



RLQ sample, we refer to Paper I. 

In order to check if some of the RQQs were detected at 
radio wavelengths, we extracted im age cutouts at the qu asar 



radio wavelengttis, we extracted im age cutouts at the qu asar 
positions from the 1.4 GHz NVSS0 ( jcondon et al. 1998| ) and 
the VLA FIRST||] survey ( [Becker, White fc Helfand 199£ 



and examined them in AlPS. We found seven quasars that 
were detected in either both or just one of the surveys and 
list the radio fluxes of these in Table 2. For the other quasars, 
we set a 3a detection as upper limit since the la noise level 
for the NVSS and FIRST are 0.5 and 0.17 mJybeam'S 
respectively. 

The quasars with radio detections are the two BQS 
quasars, three of the LBQS quasars and two of the faint 



* NRAO/VLA Sky Survey 

t Faint Images of the Radio Sky at Twenty-centimeters 
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Table 1. The quasar sample and log of observations. All images were obtained using the HiRAC equipped 
with a Ik SiTe CCD with pixel scale 0.176 arcsec, unless otherwise is specified in footnotes. In the 6th 
column, we list the FWHM of the seeing in the combined image (in the longest wavelength filter when 
there are observations in two filters). 



Quasar 




z 




Filter 


Exp. time 


Date 


Seeing 


Remarks 
















(arcsec) 




LBQS 0007- 


0003 





698 


I 


4x600 


96/07/23 


1.1 


photometric 


LBQS 0020- 


0300 





580 


V, R 


4x600 


96/07/25 


0.8 


dust 


BFSP 12344 


-00210 





753 


I 


5x300 


98/04/28 


0.7 


variable transparency 










I 


6x300^ 


00/01/03 




photometric 


BFSP 12355 


+01070 





720 


I 


9x300 


97/05/14 


0.5 


photometric 


BFSP 12364 


-0053'^ 





727 


I 


8x300 


99/05/19,23 


0.9 


variable transparency 










I 


1 X 300W 


00/01 /03 




LJ1±U LUlllL Li IL, 


BQS 1333+176* 





554 


R 


4x600 


00/01 /OS 


0.7 


photometric 


BFSP 15152 


+02441 





608 


R 


A and 


ci'7 1 r\ / 1 /I 
y / /UO/ 14 


0.6 


photometric 


BFSP 15196 


+0220 


u 


/ 4z 


T 

1 


A \y a(\r\ 

4X dUU 


riR /nv loA 
yo/u / /z4 


n Q 
U.o 


dust 


BFSP 15199 


+0247 


u 


/ Do 


-ft,, 1 


X DUU 


yo/ u / / zo 


u.o 


photometric 


BFSP 15202 


+0236 





765 


R, I 


4x600 


96/07/22 


0.6 


photometric 


BQS 1538+477 





770 


R, I 


4x600 


96/07/21 


0.6 


photometric 


BFSP 22011 


-1857 





615 


V, R 


4x600 


96/07/21 


1.0 


photometric 


LBQS 2235+ 


0054 





529 


V, R 


4x600 


96/07/21 


1.0 


photometric 


LBQS 2238+ 


0133 





714 


R, I 


4x600 


96/07/22 


0.7 


photometric 


LBQS 2239- 


0055 





680 


R, I 


4x600 


96/07/23 


0.9 


photometric 


LBQS 2245- 


0055 





801 


R, I 


5x600 


96/07/24 


1.3 


dust 


LBQS 2348+ 


0210 





504 


V, R 


4x600 


96/07/24 


1.3 


dust 


LBQS 2348+ 


0148 





749 


I 


4x600 


96/07/25 


0.8 


dust 


LBQS 2350- 


0012 





561 


V, R 


4x600 


96/07/23 


0.9 


photometric 


LBQS 2353- 


0153 





672 


R, I 


4x600 


96/07/22 


0.6 


photometric 



*Imaged using ALFOSC with a 2k Loral CCD with pixel scale 0.189 arcsec. 
timaged with the HiRAC camera, but with a 2k Loral CCD with pixel scale 0.11 arcsec. 
t Exposures taken during photometric conditions in order to calibrate the other images obtained under 
less photometric conditions. 



BFSP quasars, and we have marked these with 'i?' in Ta- 
ble 2. One of these quasars, BQS 1538+447, looks resolved 
on the FIRST map and appears to have FRI morphology. 

Can these seven (out of 20) radio-detected quasars be 
classified as RQQs? The distinction between RLQs and 
RQQs is not clear-cut, and different definitions exist. Indeed 
it seems likely that there is a continuum in radio-loudness, 
with RLQs representing the extreme of the distri bution in 
radio luminosity at a given absolute magnitude (White et 
al. 200(|). Nevertheless, a distinction is often made based 
on radio luminosity and/or the ratio of radio to optical 
flux density, R. According to the radio luminosity crite- 
rion, a quasar is radio-loud if Psghz > 10^* WHz~^sr~^, 
or PsGHz > 1-3 X lO'^^ WHz~^, wherea s acco rding to the 
criterion used by e.g. Kellermann et al. (198!:), a quasar is 
radio-loud if R <; 10. A much used definition of R is the 
ratio of the rest frame fiux at 5 GHz to the B-band flux. 



and here we choose Sramek & Weedman's (198C) definition, 
R = 55GHz/S2500, which uses the flux density at rest frame 5 
GHz and 2500 A (corresponding to B-band when z ~ 0.7). 
The usefulness of th e 7?-parameter is not clear, e.g. Gold- 
schmidt et al. (199!:) argue that R is physically meaningful 



only if the radio and optical luminosity is linearly correlated, 
and that there are various reasons that this might not be the 
case (Goldschmidt et al. and references therein). 

To calculate the _R-parameter, we assumed that both 
radio and optical fiux density scales as 5*1^ oc 1^"°, and com- 



puted rest frame flux densities at 5 GHz and 2500 A. The 
rest frame flux density at 5 GHz was found according to 



•S'sc 



5*1.4 



where we used the NVSS flux as S'i.4ghz when available, 
since this is most sensitive to extended structures. The radio 
power was also calculated, 

AgHz ~ 47rdL S'sGHz, 

where is the luminosity distance to the quasar. The opti- 
cal fiux density at rest frame 2500 A was found from the m b 
magnitudes of the quasars (determined from the images, see 
details in Section 0) using the relation 



/4400\' 



V 25007 

where logSs = -22.35 -0.4mB (|Longair 198l|). In Table 2, 



we list the rest frame radio and optical fiux densities, the 
PsGHz parameter and the radio power. 

In the above calculations, we assumed a spectral index 
of a = 0.5 for both the radio and the optical spectra, but 
for four of the quasars we estimated radio spectral indices 
by extrapolating between the NVSS/FIRST fiux densities 
and_flux densities measured at 5 GHz by Kell ermann et al. 
( tl989|) or at 8.4 GHz by Hooper et al. ( |l99^ ), denoted Or 
in Table 2. We found that BQS 1538+477 has a steep ra- 
dio spectrum with af 4 = 1.02, but the spectrum fiattens 
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Table 2. Radio fluxes at 1.4 GHz from NVSS and FIRST image cutouts. Quasars with radio detec- 
tions in NVSS or FIRST (or both) are marked with an R in the first column. For the other quasars, 
we give upper Umits that correspond to 3cr detections. Rest frame fluxes at 5 GHz and 2500 A are 
denoted Ssghz a^nd 'S'25001 a^nd their ratio is R^gHz- The radio power is denoted PsgHzi and is 
the radio spectral index between 1.4 and 5 GHz. 



Quasar 


IN V 00 


r irto i 


c .... 




S2500 


E> .... 
^SGHz 


^5GHz 






(mJy) 


(mjyj 




(mJy) 




■\A^ TT„— 1 

VV xlz 


LBQS 0007—0003 


< 1.5 


< 0.5 


<0.61 


0.50 


0.10 


< 6.1 


< 1.6x10^* 


LBQS 0020—0300 


< 1.8 


no map 


<0.76 


0.50 


0.21 


< 3.6 


< 1.4x10^* 


BFSP 12344— 0021^^^ 


4.8 


5.2 


1.92 


0.50 


0.03 


64.0 


6.1x102* 


BFSP 12355+0107 


confused 


< 0.5 


<0.20 


0.50 


0.02 


< 10.0 


< 5.8x10^^ 


BFSP 12364-0053 


< 1.5 


< 0.5 


<0.60 


0.50 


0.01 


< 60.0 


< 1.8x10^4 


BQS 1333+176^ 


37.6 


no map 


18.54 


0.320 


0.60 


30.9 


3.0x10^5 


BFSP 15152+0244 


< 1.5 


< 0.5 


<0.63 


0.50 


0.02 


< 31.5 


< 1.3x10^4 


BFSP 15196+0220 


< 1.5 


< 0.5 


<0.60 


0.50 


0.03 


< 20.0 


< 1.9x10^"' 


BFSP 15199+0247 


< 1.5 


< 0.5 


<0.60 


0.50 


0.03 


< 20.0 


< 2.0x10^4 


BFSP 15202+0236^' 


< 1.5 


1.0 


0.40 


0.50 


0.01 


40.0 


1.3x1024 


BQS 1538+447^ 


95.2 


85.4 


26.28 


1.02* 


1.00 


26.3 


8.8x10^^ 


BFSP 22011-1857 


< 1.5 


no map 


<0.62 


0.50 


0.05 


< 12.4 


< 1.3x10^4 


LBQS 2235+0054^^ 


< 1.5 


0.83 


0.94 


-0.65[| 


0.02 


47.0 


1.4x1024 


LBQS 2238+0133 


< 1.5 


no map 


<0.61 


0.50 


0.22 


< 2.8 


< 1.7xl024 


LBQS 2239-0055 


< 1.5 


< 0.5 


<0.61 


0.50 


0.12 


< 5.1 


< 1.6xl024 


LBQS 2245-0055^^ 


< 1.5 


1.3 


0.51 


0.50 


0.30 


1.7 


1.9x102-4 


LBQS 2348+02101^ 


41.8 


no map 


25.05 


0.12t 


0.17 


143.4 


3.3x1025 


LBQS 2348+0148 


< 1.5 


no map 


<0.60 


0.50 


0.08 


< 7.5 


< 1.9xl024 


LBQS 2350-0012 


< 1.5 


< 0.5 


<0.64 


0.50 


0.32 


< 2.0 


< 1.1X1024 


LBQS 2353-0153 


< 1.5 


< 0.5 


<0.61 


0.50 


0.57 


< 1.1 


< 1.5X1024 



*Spectral index found between flux at 1.4 GHz and flux measured at 5 GHz by Kellermann et al. 
(1989) 

tSpectral index between 1.4 GHz flux and flux measured at 8.4 GHz by Hooper et al. (1996) 



at high frequency, with Qg'y = —0.19 (Falcke, Sherwood & 
Patnaik 1996), so it is therefore likely that the bright core is 
flat spectrum and possibly variable, but the extended emis- 
sion is steep-spectrum. The spectrum of LBQS 2235+0054 
is somewhat unusual with af f = —0.65, i.e. the flux den- 
sity rises with increasing frequency. This could be due to 
synchrotron self-absorption at low frequencies or perhaps 
thermal absorption. 

As seen in Table 2, six of the seven radio-detected 
quasars have R > 1 0, and classify as radio-loud according to 
Kellermann et al 's ( |l989| ) criterion, but only three of these, 
the two BQS quasars and LBQS 2348+0210, have radio lu- 
minosities greater than 1.3x10^^ Wlfz^^ and qualify both 
criteria. 



Miller et al. (1993), in a study of BQS quasars, found a 
group of quasars with radio luminosities at 4.8 GHz in the 
range ~ 1.3-4.0 x 10^"* WHz"^ {Ho = 50 kms"^Mpc-\ 
qo — 0) that did not tightly follow the emission line-radio 
luminosity correlation, and termed these 'radio-intermediate 
quasars'. It was suggested that they could be Doppler- 
boosted RQQs with Lorentz factors of ~ 5. In terms of 
the i?-parameter, it is common to classify quasars with 
1 < ^?5GHz < 100 as radio-intermediate and quasars with 
^5GHz < 1 as radio-quiet. According to this, all but one of 
the seven radio-detected quasars are radio-intermediate, the 
exception being LBQS 2348+0210 with 7?5GHz = 143.4. 

Falcke et al. (1996) find that the division line at i? = 10 



spectrum, core-dominated quasars in the BQS. According 
to orientation-based Unified Schemes these are relativisti- 
cally beamed radio sources with extended lobe emission, and 
should therefore be rare and have radio flux densities higher 
than the steep-spectrum sources. Falcke et al. term these 
radio-intermediate quasars, and also suggest that they are 
relativistically boosted radio-weak quasars. They propose to 
use both the i?-parameter and the spectral index to divide 
between RLQs and RQQs, and that the division line be- 
tween radio-loud and radio-quiet is set at R = 250 for flat- 
spectrum sources and at = 25 for steep-spectrum sources, 
so that the radio-intermediate quasars have /^-parameters 
in the range 25-250. 

Blundell & Beasley ( 199^ ) also point out that if the 
criterion for radio-loudness is based on the extended radio 
emission in quasars, the radio-intermediate quasars would 
classify as radio-quiet. This would however not be the case 
for BQS 1538+447, since its extended radio flux is 53.9 mJy 
(total NVSS flux minus peak FIRST flux), corresponding to 
5.0x1025 WHz-i at 5 GHz. 

The seven quasars discussed above clearly have ra- 
dio powers much lower than typical RLQs (~ lO^^'-lO^^ 
WHz~^), and if not radio-quiet, may classify as either radio- 
intermediate or radio-weak. However, we have checked that 
our results are not affected by the presence/absence of these 
quasars. 



between RLQs and RQQs leads to inconsistencies with the 
Unified Scheme, since they find a high fraction of flat- 
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Figure 1. Absolute B magnitudes of the quasars as a function 
of redshift. Diamonds: bright BQS quasars, circles: intermediate- 
luminosity LBQS quasars, triangles: faint BFSP quasars. 



3 OBSERVATIONS AND DATA REDUCTION 

Most of the data were obtained using the High Resolution 
Adaptive Camera (HiRAC) at the 2.56-m Nordic Optical 
Telescope (NOT) during 1996 July 21-25. The HiRAC was 
equipped with a Ik SiTe CCD with a pixel scale of 0.176 arc- 
sec, giving a field of view of 3x3 arcmin. Two quasar fields 
were imaged in May 1997, also using the HiRAC, but with 
a Loral 2k chip with pixel scale 0.11 arcsec, giving a 3.7x3.7 
arcmin field of view. Images of three more quasar fields were 
obtained during observing runs at the NOT in 1998, 1999 
and 2000 using the ALFOSC (Andalucia Faint Object Spec- 
trograph and Camera) in imaging mode, equipped with a 2k 
Loral CCD with pixel scale 0.189 arcsec, and a field of view 
of 6.5x6.5 arcmin. The bulk of the data were obtained under 
photometric conditions, and the seeing FWHM is less than 
one arcsec in 15 out of 20 fields. 

Since we are attempting to detect galaxies that are 
physically associated with the quasars, the filters were cho- 
sen to give preference to early-type galaxies with strong 4000 
A breaks at the quasar redshifts. The spectral energy dis- 
tribution of early-type galaxies is characterized by a strong 
break at 4000 A with most of the energy emitted at wave- 
lengths longer than 4000 A. For an early-type galaxy at 
z > 0.67 the 4000 A break moves from R into /-band, so 
for the z > 0.67 quasars we used I-band imaging and for 
the z < 0.67 quasars we chose i?-band. Twelve of the fields 
were also imaged in two filters, either V and 7?, or R and 
I, depending on the redshift of the quasar so as to straddle 
the 4000 A break. 

Typically, the integrations were divided into four expo- 
sures of 600 s each, offsetting the telescope 10 arcsec between 
each exposure to avoid bad pixels and cosmetic defects from 
falling consistently on the same spots. The 7-band images 
from 1997-2000 (three fields) have a prominent fringing pat- 
tern due to interference of night sky lines. We therefore took 
9-11x300 s exposures of these fields in order to obtain a 
sufficient number of images to construct a fringe frame. The 
fringe frame was made by taking the median of nine /-band 
images taken close in time, and used to fiatten the back- 
ground with by dividing out the fringes. Apart from this 
fiat fielding technique to remove fringes, the data were re- 



duced as described in Paper I using standard iraf reduction 
tasks. 

We also observed standard star fields during each night 
in order to calibrate magnitudes to the standard Johnson 
photometric system. During the 1 996 an d 1997 runs, we ob- 
served fields by Christian et al. (1985), and for the other 



observing runs we selected standard stars by Landolt (1992). 

The two fields from 1998 and 1999 (BFSP 12344-0021 
and BFSP 12364—0053) were imaged as back-up targets dur- 
ing another programme that required photometric and sub- 
arcsec seeing conditions, so the transparency is somewhat 
reduced in these frames. We therefore took exposures of the 
same fields under photometric conditions in January 2000 
to enable a proper photometric calibration. 

For photometry and object detection we processed the 
images in FOCAS (Fa int O bject Classification and Analysis 
System, e.g. Valdes (198!:)), as described in Paper I. Each 
detected object was classified by using a template PSF (a 
15x15 array) formed from 4-6 selected stars in each image, 
and a record of FOCAS total magnitudes and positions of 
the detections were made by excluding objects classified as 
'multiple', 'long' or 'noise'. The detections were examined by 
eye and objects coinciding with spikes of bright stars were 
removed from the records. 

As described in Paper I, we performed completeness 
simulations in the images, and these simulations show that 
the data are complete down to 24.0 in V, 23.5 in R and 23.0 
in /, with errors of ±0.3 mag at the limits. 

The quasar fields lie at high galactic latitudes 46° < 
\b\ < 76° and the galactic reddening is E {B — V) < 0.063 in 
all fields. We corrected for galactic extinction usi ng an elec- 
tronic version of the maps by Burstein & Heiles (1982) and 
t he G alactic extinction law by Cardelli, Clayton & Mathis 
( [l989| ). 



Unless the quasar was saturated in the image, we deter- 
mined its apparent magnitude in either R or I, and there- 
after converted to B assuming an optical power-law spec- 
trum, Sv oc v'" , with a — 0.5 and zero points in B, R 
and / as given by Longair (1981). The B magnitudes were 
thereafter converted to absolute magnitudes by assuming a 
/f-correction equal to 2.5 (a — 1) log (1 + z). 

Three of the quasars were saturated in the images (BQS 
1333+176, BQS 1538-^477 and LBQS 2348-H0210) so we 
used apparent B and Bj magnitudes as given in the cat- 
alogues, ex cept f or BQS 13334-176 for which Yee, Green & 
Stockman (1986) have determined an r magnitude. The Bj 



magnitude for LBQS 2 348-1-0210 was conver ted to B by as- 
suming B = + 0.09 ([Metcalfe et al. 199l|), and the Gunn 



r magnitude of BQS 1333+176 was first converted to R by 
assuming r ^ R + 0.43 + 0.15 {B - V) ( [Kent 198E| ) and a 
_B — colour of 0.3 (i.e. _B — V of a quasar with optical 
spectral index a = 0.5), and thereafter to B as described 
above. 

Fig. |l] shows the distribution of the quasars in the z-Mb 
plane, and the Mb magnitudes are also listed in Table ^ 

3.1 Control fields 

Since we are concerned with investigating if there is an ex- 
cess of galaxies in the quasar fields, we aim to have a good 
determination of the background counts. For this purpose, 
we obtained images of random fields in the sky at approx- 
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imately the same galactic latitudes as targets in our AGN 
sample. The control fields were described in Paper I, so here 
it suffices to mention that these fields were imaged as part of 
our AGN environment survey, and therefore have the same 
depth and were obtained in exactly the same manner as the 
quasar fields. This is important in order to get a rob ust es- 
timate of the ga laxy richness in the quasar fields (Yee & 
Lopez-Cpruz 19991). 



In Paper I, we averaged the counts in the control fields 
in order to obtain a good determination of the background 
counts, and we also found that the slopes of the background 
galaxy number counts agreed w ell wi th that found by other 
investigators, e.g. Small et al. (1995). Due to the clustered 



nature of field galaxies, th e errors in the galaxy cou nts are 
expected to be ~ 1.3\/iV (fee & Lopez-Cruz 1999). How- 
ever, the variation in the counts from control field to control 
field is larger than this, presumably due to the small areas 
surveyed in each field (typically 3x3 arcmin), but also due 
to a variation in the counts with galactic latitude, especially 
in the 7-band control fields. 

This can be seen in Fig. |2| where we have plotted the 
number of galaxies brighter than the completeness limits 
from raw counts in the control fields as a function of galactic 
latitude, with error bars. The counts in the quasar 

fields within a 0.5 Mpc radius centered on the quasar are 
also plotted, and we have included the RLQ fields from Pa- 
per I. The error bars on the counts in the quasar fields are 
intentionally left out for clarity. It is seen that six of the 
eight /-band control fields lie at galactic latitudes between 
45-60°, where also the bulk of the quasar fields are found. 
The remaining two lie at |&| ~ 30°, and were obtained in 
order to get a handle on the counts here since two fields in 
our radio galaxy sample (to be published in a future paper) 
lie at lower latitudes. 

The counts in the two \b\ ~ 30° fields have a mean 
of 16.94±2.3 arcmin"^, whereas the fields between 45 and 
60° have a mean of 11.53±1.6 arcmin"'^, so the background 
counts appear to increase towards lower latitudes. Presum- 
ably this is due to faint, red stars as the background counts 
in R do not seem to vary as much with latitude. It is seen 
that the bulk of the quasar and control fields lie at latitudes 
between 40 and 60°. Nevertheless, in addition to using the 
average background counts, as we did in Paper I, we also 
here make a correction for the variation in the counts with 
latitude. 

In /-band, a least square fit to the eight data points 
weighted with the error bars gives a relation for the back- 
ground counts of N (/) = -0.32|fe| -I- 26.6. This is shown 
by the solid line in Fig. ^, and we call this Model 1. A 
straight line through the mean values of the two low latitude 
and the six high latitude fields gives a less steep relation of 
(/) = -0.25|6| -I- 24.5, as shown by the dotted line. This 
relation is named Model 2. 

As mentioned above, the /?-band background counts 
show a weaker dependence on galactic latitude. The cor- 
responding Model 1 for the R counts is TV (R) = — 0.15|&| -f- 
18.6, and by drawing a line through the mean of the seven 
fields at 40-60° and the mean of the two |fe| « 30° fields, 
we get that Model 2 for the /?-band counts is A'^ (/?) = 
— 0.17|6| -I- 19.4 galaxies per arcmin'^. Outside the ranges of 
the sloped lines in Fig. H, we have used the constant values 
corresponding to the endpoints of the relations. 



We have chosen to let the background counts follow 
these relations when correcting for background galaxies in 
the quasar fields, but we have also used the average of all 
control fields, named Model 3. Details about these models 
can be found in Table ^. We have also included the V^-band 
control fields in this table since these were used to find the 
excess in five quasar fields that were also imaged in R. Since 
there are only two V-band control fields (at \b\ = 57?4 and 
46? 7), we have used as Model 1, a relation that has the same 
slope as the R counts, but is normalized to the mean V'-band 
counts. 



4 ANALYSIS 

To quantify the excess counts in the quasar fields we have 
used the amplitude, /3gq, of the spatial galaxy-quasar cross- 
correlation function, ^ (r) — B^qT where 7 = 1.77. The 
amplitude is evaluated at a fixed radius of 0.5 Mpc at the 
quasar redshift, corresponding to « 1 arcmin at z = 0.7 and 



has units of Mpc ' . Longair & Seldner (1979) showed how 



it can be used to estimate the clustering about any point in 
the Universe using galaxy counts, and it has subsequently 
been widely used to quantify galaxy environment around 
diflferent types of AGN (e.g. Yee & Green 1987; Yates, Miller 
& Peacock 1989; Ellingson et al. 1991a; Smith, O'Dea & 
Baum 1995; Wurtz et al. 1997; DeRobertis, Yee & Hayhoe 
1998; McLure & Dunlop 2000). In Paper I we used it to 
quantify the environments around RLQs, and here we apply 
the same analysis to the RQQs. We therefore refer to Paper I 
for a more detailed description of the analysis, and give only 
a brief summary here (see also Longair & Seldner 1979; Yee 
& Green 1987). 

First, the amplitude, ^gq, of the angular cross- 
correlation function is found. This is easily obtainable from 
the data since it is directly proportional to the number of 
excess galaxies within the 0.5 Mpc radius centered on the 
quasar. Thereafter, /3gq is obtained from the following ex- 
pression 

A/'gAgq 



$ (miim, z) L, 

where A/'g is the average surface density of background galax- 
ies, de is the angular diameter distance to the quasar and 
I-, = 3.78 is an integration constant. 

In the denominator, $ (miim, z) is the luminosity func- 
tion integrated down to the completeness limit of the data, 
thus giving the number of galaxies per unit volume at 
the quasar redshift seen in the fiux-limited data. We have 
used a Schechter function to evaluate the integral, with 
characteristic magnitude M^g (B ) — —20.83 and slope 
Q — —0.89 based on Lilly et al's (1995) estimation of the 
CFRS (Canada-France Redshift Survey) luminosity function 
in the 0.5 < z < 0.75 redshift interval. Since the integrated 
Schechter function is proportional to the characteristic den- 
sity, (p* , the /3gq estimator becomes inversely proportional 
t o it, and therefore very sensitiv e to it s value. Yee & Green 
(1987) and Yee & Lopez-Cruz (1999) have advocated the 



importance of choosing (j)* in such a way that it is consis- 
tent with the data. In Paper I, we explain how we used the 
CFRS luminosity function to construct galaxy counts that 
were fitted to the observed counts in the control fields. The 
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Table 3. Various relations used for the galaxy background counts. In Model 1 and 2 the counts 
are dependent on the absolute value of the galactic latitude |b|, whereas Model 3 is just the average 
counts. The error in the average given in the last column is the standard deviation in the mean, 
and the other errors were calculated as l.SvW- 

Model 1 Model 2 Model 3 



Filter 


Range (°) 


TV (arcmin 


Range (°) 


N (arcmin ^) 


N 


(arcmin 


/ 


30.0< |fe| <57.1 


-0.32|6|+26.6 


30.1< |b| <51.6 


-0.25|fe|+24.5 


12 


36±4.38 


/ 


|b| <30.0 


17.07±5.37 


\b\ <30.1 


16.94±5.35 


12 


36±4.38 


I 


\b\ >57.1 


8.50±3.79 


\b\ >51.6 


11.53±4.41 


12 


36±4.38 


R 


29. 2< |6| <57.2 


-0.15|6|+18.6 


29. 7< |b| <48.5 


-0.17|6|+19.4 


12 


51±2.69 


R 


|b| <29.2 


14.07±4.88 


\b\ <29.7 


14.45±4.94 


12 


51±2.69 


R 


\b\ >57.2 


9.74±4.06 


\b\ >48.5 


11.30±4.37 


12 


51±2.69 


V 


30.0< |b| <57.4 


-0.15|6|+18.23 






10 


19±6.88 


V 


|b| <30.0 


13.60±4.79 






10 


19±6.88 


V 


\b\ >57.4 


9.40±4.0 






10 


19±6.88 



— band fields R — band fields 




20 30 40 50 60 70 80 90 30 40 50 60 70 80 90 
bl (degrees) Ibl (degrees) 



Figure 2. The filled circles show the counts of galaxies in the /- and ij-band control fields as a function of galactic latitude (absolute 
value). The number of galaxies in the quasar fields are also plotted using asterisks for the RQQs and triangles for the RLQs. The solid 
and dotted lines show the relations for the background counts as a function of galactic latitude, corresponding to Model 1 and Model 2, 
respectively (see text and Table |^ for details). 



fitting gave (pl = 0.0065, (j}*^ = 0.0072 and 01 = 0.0052 
Mpc~^, and we use these numbers in our computation of 

$(miim,2). 



Yee & Lopez-Cruz (1999) have examined the robust- 



ness of -Bgc (index 'gc' means galaxy-cluster centre) on a 
sample of 47 low-z Abell clusters, and find that it is a ro- 
bust estimator of galaxy richness in clusters. They apply a 
number of tests, and find that Bgc is fairly robust provided 
the normalization of the luminosity function reproduces the 
observed galaxy counts. Under this condition, varying e.g. 
M* causes variations in Bgc of ~ 10 per cent. Also counting 
dovirn to different absolute magnitude limits produces sys- 
tematic variations of only a few per cent, but the best and 
most reliable results are obtained if one counts down to 1-2 
magnitudes fainter than M* . In the J-band images we have 
reached approximately 1.7 mag past M* , and in the _R-band 
data, about 2 mag (see last column of Table 0). 



5 GALAXY EXCESS IN QUASAR FIELDS 

We counted the number of galaxies within the 0.5 Mpc ra- 
dius in the quasar fields and averaged the counts for the 
R and /-band data, i.e. for the seven z < 0.67 and the 
13 z > 0.67 quasar fields. The average counts are shown 
in Fig. ^ where we also have plotted the average galaxy 
counts in the control images. The error bars of the back- 
ground galaxy count were calculated as 1.3y/N. There is 
a clear excess of faint galaxies at 7 > 20, but the galaxy 
counts in the 7?-band fields show no average excess above 
the background. 

For each quasar field, the net excess of galaxies within 
the 0.5 Mpc radius, A'^net, was calculated by subtracting the 
background counts according to Model 1-Model 3 from the 
counts in the quasar fields. Thereafter, the amplitudes 
and Bgq were obtained. The results are shown in Table m 
where we, for each model of the background counts, show 
the net excess, Anct, the significance, a, of the excess, and 
the clustering amplitude, Bgq. 

Since the background counts were scaled to the area of 
the 0.5 Mpc circle, TVnct is not an integer. The significance. 
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Figure 3. The open circles sliow the average galaxy counts within 
the 0.5 Mpc circle in the quasar fields, and the filled circles show 
the counts in the control fields. The vertical dotted lines mark 
the completeness limits, / = 23.0 and R = 23.5. 



a, was calculated as Nuet/l-3\ 



and the uncertainty in 



A^gq _ (iVnct + 1.3^ iVb) 
-Bgq A'nct 



1/2 



( [Yee fc Lopez-Cruz 1999| ). 

In the last column of Table ^, we show approximately 
how many magnitudes fainter than M* an Sa galaxy with 
apparent magnitude m = is. This demonstrates that 

we are close to the recommended M* + 2 in all the fields so 
that a reliable estimate of Bgq is obtained. By transforming 
M*AR (B) = -20.83 to R and I, assuming Bab = B - 0.17 



(Oke 1972) and using colours of a zero redshift Sa galaxy 
r nodel, B - V = 0.74, V - R = 0.68 a nd R - I = 0.57 
(Guiderdoni & Rocca-Volmerange 1988), we found that 



M*n = -22.08 and Mj = -22.65. 

The mean clustering amplitude using Model 1 for 
the background counts is 336±77 Mpc^'^*^, whereas using 
Model 2 and Model 3 gives lower mean amplitudes of 212±74 
and 210±82 Mpc^ '^'^. The correction in the background 
counts to account for galactic latitude is therefore more dras- 
tic in Model 1, giving a higher mean Bgq, but the three mean 
values are nevertheless within the errors of each other. (The 
errors we give for the mean Bgq's are standard deviations 



Table ^ are consistent with each other, and their variation 
should give some indication of the real errors that are in- 
volved in this analysis. 

The mean clustering amplitude in the quasar fields, re- 
gardless of which model is used for the background counts, 
is higher than the local galaxy-galaxy correlation amplitude 
of 60-70 Mpc^ ''^ (Davis & Peebles 1983; Loveday et al. 
1995; Guzzo et al. 1997), so on average, we find that the 
RQQs are not located in field-like environments, but rather 
seem to prefer groups or poorer clusters. 

Using a test, we compared the distribution of the 
significances of the excess to that of a normal distribution 
(which is what we expect if there are no excess galaxies 
around the quasars), and find a probability of 9x10^'^" 
(Model 1) that the distribution is equal to that of a nor- 
mal distribution. For Model 2 and 3, we get probabilities of 
9x10"^ and 1x10"^°. Excluding the seven radio-detected 
quasars, we get probabilities (for Model 1 to 3) of l.Ox 10"^", 
1.1x10^^ and 4.4xl0~^ that the distribution of cr's is equal 
to a normal distribution. Furthermore, excluding the two 
BQS quasars, the probabilities become 5x10"^", 5xlO~* 
and 7x10"''. These very low probabilities agree with the vi- 
sual impression from Fig. ^, where the Bgq's of both RLQs 
and RQQs are plotted as a function of redshift and Mb ■ The 
low probabilities are a strong result, and can imply either 
an average density enhancement or a variance in the envi- 
ronments greater than expected, and most likely, both are 
present in the data. 

This adds more weight to the result that RQQs are 
on average not located in field-like environments. However, 
the environmental richness varies between individual fields; 
some fields have no significant excess, whereas other appear 
to be very rich in galaxies. There are five fields with a > 2a 
excess irrespective of how the background is calculated. The 
richest field among these is the field around BQS 1538-1-447 
with a clustering amplitude in the range 1100-1200 Mpc^ '^'^ 
and a galaxy excess of ~ 40-50. In Paper I we argued that 
an amplitude of ~ 740 Mpc^'^^ corresponds roughly to Abell 
richness class ^ 1, so the cluster candidate around BQS 
1538+447 may even qualify as an Abell class 2 cluster. Four 
other fields have galaxy excesses of 20-30 galaxies and am- 
plitudes ^ 500 Mpc^ '^'^ and are probable Abell class 1 clus- 
ters. BFSP 151994-0247 and 15202+0236 have amplitudes in 
the range 700-800 Mpc^'^'', whereas LBQS 2239-0055 and 
2353—0153 have somewhat smaller amplitudes of 400-500 
Mpc^-^". 

As Fig. ^ shows, the _R-band fields have a galaxy surface 
density close to that in the control fields, and the mean clus- 
tering amplitude for these, 71±34 Mpc^ '^^, is consistent with 
that of galaxies in the local universe (23±50 and — 43±50 
Mpc^ ''^ for Model 2 and 3). For the 13 7-band fields the 
mean clustering amplitude is 478±96 Mpc^ '"'' (314+101 and 
346+106 Mpc^ " for Model 2 and 3, respectively). It might 
therefore seem that there is an apparent trend with redshift 
since the J?-band fields represent the z < 0.67 quasars, but 
the current sample of 20 objects is too small and span a too 
narrow redshift range to test this properly. 

In principle, we are interested in the full distribution 
function for the clustering amplitude, e.g. when the data are 
compared with models. The distribution will be positively 
skewed with a long tail to the right since we do not observe 
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Table 4. The results of the analysis of galaxy counts in the RQQ fields. For each of the three models for the background counts 
(see Table the net excess, A^nct, the significance of the excess, cr, and the clustering amplitude, -Bgq, are listed. The absolute 
B magnitude of the quasar is Mg, and AM in the last column shows how many magnitudes into the luminosity function the 
completeness limit extends, i.e. M* + AM. All numbers are based on galaxy counts in either R or I depending on whether the 
quasar redshift is < 0.67 or > 0.67. 
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fields with negative Bgq's comparable in absolute value to 
the most positive Bgq values. For a positively skewed distri- 
bution, the median will always be smaller than the mean, 
and for our sample the median _Bgq is 131 Mpc^'^^ (Model 3). 
A model for formation of RQQs should be able to reproduce 
both the mean amplitude and the intrinsic scatter in Bgq. 
Thus the intrinsic scatter is equally important as the mean 
amplitude. The scatter (standard deviation) of the Bgq dis- 
tribution is asgq = 365 ± 58 Mpc^'^^ (Model 3), and this 
includes both intrinsic scatter and measuring uncertainty. 
We have used a measuring uncertainty of ABgq ~ 185 ± 48 
Mpc^'^^, and estimated the intrinsic scatter, crint, assuming 
that o-f„t = + AB^q, giving aint = 315 ± 73 Mpc^'^^. 

As expected, the intrinsic scatter is large, but difficult to 
estimate due to the large measuring uncertainty in Bgq. In 
Table ^ we list values of (Tint, mean and median Bgq's that we 
get by using the different models for the background galaxy 
counts. 

Since 12 of the fields were imaged in two filters, we 
also computed Bgq in the fields taken in the shortest wave- 
length filter, and list the results in Table |^. The mean am- 
plitudes are 220±92 and 120±93 Mpc^'^^ for Model 1 and 
Model 3, respectively (Model 2 gives values intermediate be- 
tween these two). For comparison, the mean amplitudes for 
the same 12 fields obtained in the reddest filters are 367±120 
and 276±125 Mpc^ ". So, the mean amplitudes are lower, 
but still within the errors of each other, and on average, the 
excess is detected in both filters. There are two fields which 
have significantly higher Bgq in / than in R, and these are 



BQS 1538+447 and LBQS 2353-0153, with Bgq approxi- 
mately a factor of two higher in 7 than in i?, suggestive of a 
fair amount of red galaxies in these fields. This is what we 
expect if there are many red, early-type galaxies with strong 
4000 A breaks at the quasar redshift. 

Although the overall trend is that Bgq is larger in the 
reddest filters, there are 3-4 fields where the opposite is seen. 
This can, apart from in one case, probably be attributed to 
measurement errors since these are fields with no significant 
excess. The exception is BFSP 22011-1857, with a ~ 2o- 
excess in 1/, but with no excess in 7?. The clustering am- 
plitude deduced from the 7? image lies in the range 50-150, 
whereas it is « 470 Mpc^'^^ in F-band. This field was care- 
fully checked in order to verify that there are more galaxies 
in the V image than in the 7? image. There are two possible 
explanations for this, one is that there is a foreground con- 
centration of bluer galaxies toward this quasar, the other is 
that it might be a cluster with a blue galaxy population at 
the quasar redshift. 



5.1 Properties of the excess galaxies 

Elliptical and early-type spiral galaxies in clusters are known 
to form a tight sequence in the color-magnitude diagram, 
the so-called 'red sequence'. It is the strongly evolved 4000 
A breaks in these galaxies that make their colour red, since 
metals in the cool atmospheres of giant stars absorb much 
of the radiation blueward of 4000 A. 

In the conventional models, elliptical galaxies are be- 
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Table 5. The results of the analysis of the quasar fields that were imaged in two filters, 
here in the shortest wavelength filter, i.e. in V for the z < 0.67 quasars and in R for the 
z > 0.67 quasars. The completeness limit in V is 24.0. 
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lieved to form during a monolithic collapse at 2 > 2 and 
evolve passively after an initial starburst (e.g. Larson 1975; 
Arimoto & Yoshii 1987). An alternative model that has re- 
ceived much attention lately is the hierarchical formation 
scenario, where elliptical galaxies are formed by mergers of 
disk systems (Toomre & Toomre 1972). It has been shown 
that this scenario is also consistent with the presence of 



the red sequence in galaxy clusters ( KaufTmann & Chariot 
199f~| 

Here, we have picked the four richest quasar fields that 
were imaged in two filters, and examined the colours of 
the galaxies. The four fields are BFSP 15199+0247, BFSP 
15202+0236, BQS 1538+447 and LBQS 2353-0153, with 
Bgq k, 500 Mpc^ ''^ and a ^ 3cr excess. The mean redshift 
for these quasars is (z) = 0.74 ± 0.02. Fig. |^ shows a colour- 
magnitude diagram of galaxies in these fields, where we have 
excluded stars at / < 21 (as explained later in this section). 

There is a hint of a red sequence ati? — /~ 1.5-1.7 
in this figure, i.e. tentative evidence that these fields likely 
contain galaxy clusters at z ~ 0.7-0.8. Very few galaxy clus- 
ters are know at z ^ 0.7, but judging from other work, the 
red sequence is expected to occur at 1.3 ^ R — I ^ 2.0. 
(Luppino & Kaiser 1997; Clowe et al. 1998, Lubin et al. 
2000) . Also spectral synthesis models predict that this is the 
colour expected for early-type galaxies at redshifts 0.7-0.8 
(e.g. Fukugita, Shimasaku & Ichikawa 1995). 

The red sequence is known to steepen toward brighter 
magnitudes, and this can possibly also be seen in Fig. ^ 
The steepening toward brighter magnitudes is believed to 
be caused by the mass-metallicity relation in the galaxies, 
the most luminous and massive galaxies having a higher 
metal content and therefore also stronger 4000 A breaks. 
The slope of the red sequence in clusters is also expected to 
flatten toward higher redshifts. As the formation epoch of a 
galaxy is approached, the stellar populations in the massive 
ellipticals become younger and bluer relative to the low-mass 
systems, and by measuring the slope of the red sequence as 
a function of redshift, one can constrain the formation epoch 
of the galaxies. Gladders et al. (199J) have studied clusters 
in the redshift range < z < 0.75, and find that there is little 
evidence for a change in the slope prior to z = 0.75, thereby 
concluding that the formation epoch occurs at z > 2. For 
two clusters at z = 0.756, Gladders et al. find a slope of 
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Figure 4. Colour-magnitude diagram of galaxies in the four rich- 
est quasar fields, showing a hint of a red sequence at R— I ^ 1.5- 
1.7. Possible stars at 7 < 21 are excluded. The dashed line 
shows a least square fit to the galaxies at 19 < / < 21.5 and 
1.3 < i? - / < 2.0 with a slope of -0.07+0.03. The soUd line 
marks the completeness limit in R. 



— 0.0566±0.0066, and a least square fit to the galaxies at 
19 < / < 21.5 and 1.3 < R- I < 2.0 in Fig. | gives a slope 
of « — 0.07±0.03, consistent with Gladders et al. 
Ferreras 



Silk (2000) also show that the red sequence 
in clusters can by used to discriminate between the two for- 
mation scenarios, monolithic vs hierarchical, particularly at 
the bright end. The two formation scenarios predict signif- 
icantly different bright-end slopes at high redshifts, and a 
number of z > 0.7-0.8 clusters must be observed to test 
this. Since few clusters are known at these redshifts, AGN- 
selected clusters are therefore useful targets for such studies. 

Can the tentative red sequence at i? — / ~ 1.5-1.7 be 
caused by stars? This is probably not the case, since we do 
not expect many stars in these high latitude fields (these 
four quasar fields lie at latitudes 46-61°). Nevertheless, the 
7-band counts in the control fields seem to indicate that red 
stars are present in the fields at some level. To investigate 
this further, we made an attempt to identify stars in these 
images by their half-light radii, as shown in Fig. ^. 

This was done using the imcat software by N. Kaiser, 
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developed to measure the distortion of faint galaxies induced 
by gravitational potentials of galaxy clusters. The software 
automatically measures the half-light radii of galaxies in the 
images, a feature that is not included in FOCAS. However, 
IMCAT uses a different detection algorithm than FOCAS. In 
FOCAS, objects are detected by locating connected regions 
above a certain threshold, whereas imcat smooths the im- 
ages with gaussian filters using a range of different radii (see 
Kaiser, Squires & Broadhurst (1995) for details). 

In the upper plot in Fig. ^, the stellar locus is seen 
along a half-light radius of ~ 0.4 arcsec, corresponding to a 
seeing FWHM of ~ 0.8 arcsec, and stars can clearly be sep- 
arated from galaxies down to 7 « 21. We selected objects 
with half-light radii between the two horizontal, solid lines 
(corresponding to half-light radii of 1.9 and 2.5 pix, or 0.33 
and 0.44 arcsec) and plotted them in a colour-magnitude 
diagram as shown in the lower plot. The colour-magnitude 
distribution of stars is expected to be bimodal, with disk 



stars at the blue end and spheroid stars at the red end (Bah- 
call & Soneira 1981). The bimodal distribution is not clearly 
seen here, presumably due to the small number of stars, but 
there seems to be a gap at 0.8 < R — I < 1.3. Bahcall & 
Soneira find that the stars in the Galactic spheroid have 
typical colours of i? — / « 1.6-2.2, and that the bluer disk 
population has R — I ~ 0.1-0.8, which appears to agree 
roughly with what we see here. 

At 19 < / < 23 and at colours l.S < R - I < 2.0, 
where we suspect the red sequence to be, 18 out of 90 ob- 
jects have half-light radii corresponding to the stellar locus 
in Fig. ^ i.e. four-five stars per field, so the red sequence 
cannot be caused exclusively by stars. Since imcat appears 
to have done a very good job at identifying the stars, we 
omitted the objects on the stellar locus at / < 21 in the 
colour- magnitude diagram in Fig. ^ The R — I colours were 
evaluated using aperture magnitudes measured by IMCAT in 
a 1.5 arcsec aperture, corresponding to approximately twice 
the FWHM of the seeing disk. 

We also examined the radial distribution of galaxies in 
the quasar fields as shown in Fig. ^ Here, we have taken 
the average of all fields, and counted galaxies in annuli of 
100 kpc, and also corrected for missing areas due to large, 
extended foreground galaxies or stars. As seen from this fig- 
ure, the distribution of galaxies seems to fall off outwards, 
but it is not completely clear whether the quasars lie at the 
centre of the distributions. 



6 COMPARISON WITH OTHER STUDIES 

In this section we compare our result with those obtained for 
other samples of RQQs. A compilation of results from the 
literature is shown in Fig. ^, where we have plotted mean 
clustering amplitudes as a function of median redshift and 
absolute quasar B magnitude. 

The lowest redshift point in Fig. ^ is the mean ampli- 
tude found for 33 CfA Seyfe rt galaxies (both type 1 and 



2) by DeRobertis et al. (1995), Bgq = 37.4 ± 12.9 Mp 



DeRobertis et al. compute mean clustering amplitudes sep- 
arately for Seyfert 1 and Seyfert 2 galaxies, but these are 
consistent with each other in the case where the counting ra- 
dius is 0.5 Mpc. Furthermore, Seyfert 2 galaxies are believed 
to contain hidden Seyfert 1 nuclei (e.g. Antonucci 1993), so 
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Figure 5. The upper plot shows the half-light radius of the ob- 
jects in the four richest quasar fields as a function of magnitude. 
The stellar locus is clearly seen as the objects between the two 
horizontal lines (i.e. at seeing FWHM 0.8 arcsec). The lower 
plot shows a colour-magnitude diagram of the objects on the stel- 
lar locus. 



we have chosen to use the mean for all Seyferts. Seyfert 1 
galaxies are the low-luminosity counterparts of RQQs, so we 
have assumed a mean nuclear luminosity of Mb = —21 for 
this sample. The median redshift is 0.0192 (from table 1, 
DeRobertis, Hayhoe fc Y ee (1998)). 



Yee & Green (1987) found clustering amplitudes for 
BQS quasars in three redshift ranges, 0.00-0.15 (17 quasars), 
0.15-0.30 (16) and 0.30-0.50 (7). They list the ratio of the 
mean amplitudes to the amplitude of the galaxy-galaxy 
correlatio n function, Bfra/B^^. Assuming that Bgg = 67.5 
Mpc^ '^^ ( [Davis fc Peebles 1983| ), we find that Bgq is 113±36, 
129±45 and 73±45 Mpc^ " in the three redshift ranges, re- 
spectively. To find the median absolute magnitude in the 
three sub-samples, we converted the appare nt r m agnitudes 
listed by Green & Yee ( |l984| ) and Yee et al. ( |l986| ) to R, and 
thereafter to B as described in Section |^. We found median 
absolute B magnitu des of -23.4, -24.24 and -25.83. 

Smith et al. (1995) report a mean amplitude of 



Bgq/Bgg = l.Olg'4 (2"" errors) for their sample of 169 
z < 0.3 X-ray selected quasars with median redshift 0.15. 
The quasar My magnitudes in this sample span a wide 
range, from ~ —27 to ~ —20, but since there is no sig- 
nificant difference between the amplitudes of the angular 
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Figure 6. Average radial surface density profiles for the quasar 
fields. The quasars are located in the inner bins, and the error bars 
show the standard deviation in the mean. The dotted horizontal 
lines correspond to the average background count. 



correlation function for their high- and low luminosity sub- 
samples, we have used the peak of the magnitude distri- 
bution. My = —22.2, converted to Mb = —21.9 assuming 
B-V = 0.3. 

Fisher et al. ( [L996| ) find an average Bgq of 72±20 
[h~^Mpc]^'^'^ for a sample of 14 optically bright RQQs, 
with a median redshift of 0.24, and a median AIv of —23.5 
{h = 1). Assuming B — V = 0.3 and h — 0.5, we transformed 
to Bgq = 246 ± 66 Mpc^-^^ and Mb = -24.71. The RQQs 
i n Fis her et al's sample are included in McLure & Dunlop's 
(2000) sample, along with seven other sources. McLure & 



Dunlop find a mean Bgq of 326±94 Mpc^- ' ' and their sample 
has a median redshift of 0.165, and a median My = —24.4:, 
which was converted to Mb = —24.1. 

At median z = 0.4, we find the 32 RQQs in Elling- 



son et al's (1991a) sample. The mean amplitude is 16±45 
Mpc^'^^, and we transformed the appar ent qua sar r magni- 
tudes (given in table 2, EUingson et al. (1991b)) to absolute 
B magnitudes as described in Section ^ and found a median 
Mb = -22.78. 

The 82 X-ray selected quasars studied by Smith et al. 



(2000) cover the redshift range 0.3-0.7, and in their figure 8, 



-Bgq is plotted as a function of Mv — 5 log ft. By reading 
pairs of Bgq and AIv — 5 log h off the axis of this figure, 
we estimate that the z < 0.5 and the z > 0.5 sub-samples 



who find a mean Bgq = 3 ± 40 
-23.75, and we have assumed 



have median My — 5 log ft ~ —21.8 and ~ —22.4, respec- 
tively. These were converted to Mb = —23.00 and —23.61 
assuming ft = 0.5 and B — V = 0.3. We found the mean am- 
plitudes in these two redshift intervals to be approximately 
6.13±2.95 and 4.52±7.88 [h-^Mpc]^ ". (We estimate that 
we can read off My — 5 log ft and Bgq with an accuracy of 
±0.1 and ±5 [h~^Mpc^'^'^], respectively.) Smith et al. use a 
definition of Bgq which assumes that the clustering is stable 
in co-movin g coo rdinates, whereas the definition by Longair 
& Seldner ( 1979 ) that we employ, assumes stable cluster- 
ing in proper coordinates. We therefore transformed Smith 
et al's amplitudes to our convention by multiplying with 
(1 + z)^~^' , assuming that 7 = 1.77, and that z = 0.4 and 
z — 0.6 for the two sub-samples. Thereafter, we transformed 
to the h = 0.5 case, and got mean amplitudes of 31.6±15.2 
and 27.5±47.9 Mpc^ '^''. 

The next point along the redshift axis represents the 
results reported by Boyle, Shanks & Yee ( 198^ ) for eight 
BFSP quasars in the redshift range 0.55 < z < 0.7, with 
mean Bgq/B^^ = 2.3 ± 1.3. We have assumed Mb = —23.5 
and a mean redshift of 0.625 for these quasars, since no other 
information is given. 

Our own sample is plotted with a mean Bgq = 210 ± 82 
(using Model 3), and a median redshift and Mb of 0.714 
and -24.40. 

Finally, the high-z point corresponds to the work by 
Boyle & Couch ( |1993| ) 
Mpc^ ". The median Mb is 
a mean redshift of 1.2. 

There are several aspects about the quasar phe- 
nomenon, such as formation, fuelling and evolution, that we 
can learn about by looking at their environments as a func- 
tion of redshift and AGN luminosity. But this has been diffi- 
cult to do, both because there appears to be a wide range in 
richnesses within samples, and also because it requires good 
coverage of the redshift-luminosity plane. A figure like Fig. ^ 
might help, since looking at the mean Bgq values over a wide 
range in redshift and luminosity is more meaningful than re- 
lying on amplitudes found in single fields. The drawback is 
of cause that Fig. ^ includes results from difi'erent studies 
with differing selection criteria, observations, analysis etc., 
and therefore contains uncontrollable selection effects and 
biasses. We are therefore careful with interpretations, and 
merely look at it as a useful way of comparing different sur- 
veys to see how they agree or disagree with each other as we 
move along the redshift and the luminosity axis. 

Approximately half of the data points are significantly 
higher than the galaxy-galaxy clustering amplitude of Bgg — 
67.5 Mpc^'^'^, whereas the other half is consistent with, or 
less than, Bgg. The samples from the large cross-correlation 
studies (Smith et al. 1995; 2000; Boyle & Couch 1993) be- 
long to the latter, along with EUingson et al's sample, Yee 
& Green's high-z sample and the Seyfert sample of DeR- 
obertis et al. The Seyfert sample is perhaps not so difficult 
to explain if we consider that spiral galaxies prefer field- 
like environments. If we further assume that low- luminosity 
quasars are mostly hosted by spiral galaxies, this might also 
explain EUingson et al's result, but it does not explain why 
the mean amplitude of Yee & Green's bright, high-z sample 
is so low. 

It is however noteworthy that the datapoints which 
agree with a galaxy environment richer than that of field 
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galaxies are clustered around Mb ~ —24. The exceptions 
are Boyle & Couch's sample and Smith et al's (2000) 0.5 < 
z < 0.7 sample, which have Bgq < Sgg- This may lead us to 
think that in these two studies, perhaps the galaxy excess 
(if present) was not detected. Their observations were done 
with filters that targets emission on the short-wavelength 
side of the 4000 A break at the quasar redshift, and we 
found in Section ^ that the clustering signal in fields with a 
clear excess is weaker shortward of 4000 A. Boyle & Couch 
used R and a wide VR passband (a combination of John- 
son V and Kron-Cousins R), and a,t z ~ 1.2 both these 
filters target emission shortward of the 4000 A break, and 
the sensitivity to early-type galaxies might not be partic- 
ularly good. Also, the limiting magnitude, R — 23, might 
be somewhat too shallow in order to reach sufficiently faint 
i nto t he luminosity function at these redshifts. Smith et al. 
(2000) used V-band, also corresponding to emission short- 
ward of the 4000 A break, at least at z > 0.4 where the bulk 
of their quasars are. 

Most of the clustering studies similar to our (in terms of 
method and analysis) seem to agree fairly well that RQQ en- 
vironments are different from field-like environments, both 
at high and at low redshifts. Yee & Green, Fisher et al. 
and McLure & Dunlop, looking at bright quasars at low 
redshifts, get clustering amplitudes consistent with ours, al- 
though Yee & Creen's is a bit lower than the other two. 
Boyle et al's mean amplitude for BFSP quasars at 0.55 < 
z < 0.7 also agree with ours at z « 0.7 and Fisher et al's 
at 2 ~ 0.2, so this raises the possibility that in fact there is 
little change in RQQ environments from z ~ 0.2 to 2 ~ 0.7- 
0.8. This is in agreement with recent studies that also find 
no changes in RLQ environments with redshift (Paper I; 
McLure & Dunlop 2000). 

In light of the above, and considering the similarities 
in method and analysis, the low amplitude from EUingson 
et al's work is difficult to explain. Since our quasar sam- 
ple cover a wide range in B luminosity, there are several 
quasars at similar luminosities to EUingson et al's quasars. 
The six quasars at Mb > —23.5 in our sample have a 
mean Bgq = 237 ± 78 (Model 3), which is very different 
from EUingson et al's 16±45 Mpc^'^^. The remaining 14 
quasars at Mb < —23.5 have a mean amplitude of 196±110 
Mpc^ " (Model 3), so we find no difference between high- 
and low- luminosity sub-samples. Since the six quasars at 
Mb > —23.5 have a median redshift of 0.720, i.e. higher 
than EUingson et al's z — 0.4, it might seem as though the 
environmental richness change with redshift for the lower 
luminosity quasars, but not for the high luminosity ones. 
This is perhaps somewhat odd since we have no reason to 
believe that our amplitudes have a systematic error relative 
to EUingson et al's. 

In Paper I, we found that the -Bgq's for our RLQ fields 
lay in the same range as those found for RLQs by EUingson 
et al., and the mean amplitudes were consistent with each 
other. Fur therm ore, for the RLQ MRC 0405-123, Yee & 
EUingson (1993) found Bgq = 905 ± 277, whereas we found 



579±179 Mpcr^ , i.e. slightly lower, but still in agreement 
with each other. So it seems unlikely that our amplitudes are 
systematically different from EUingson et al's. One possible 
explanation for the discrepancy may be the inhomogeneity 
of EUingson et al's sample. The ir sample was chosen from 
the Hewitt-Burbidge catalogue ( Hewitt fc Burbidge 1987 ), 



which is an incomplete compilation of quasars, but neverthe- 
less the best catalogue available at that time. There could 
be some subtle bias in this sample that causes the highest 
_Bgq fields to be excluded. E.g. if we clip out the five high- 
est -Bgq's in our sample, we get a similar result as EUingson 
et al., most RQQs having Bgq less than 300 Mpc^'^'^ and a 
mean consistent with their. 

We have seen that, whereas some surveys find that 
RQQs live in field-like environments, others agree with our 
survey finding a richer than average galaxy density around 
RQQs. But the question whether the disagreement is due to 
selection effects caused by redshift and/or luminosity, ap- 
pears difficult to answer. 



7 A DIRECT COMPARISON WITH RLQ 
ENVIRONMENTS 

In Paper I, we investigated the galaxy environments around 
21 steep-spectrum RLQs. The current sample of RQQs 
matches the RLQ sample in redshift and B magnitude, and 
both samples were observed with the same telescope and 
analysed in the same manner. This offer us a direct and 
internally consistent way of comparing the galaxy environ- 
ments around RLQs and RQQs at the same epoch and with 
similar AGN luminosities. 

In Fig. ^, we show Bgq as function of z and AIb for 
both samples (using Model 3 for the background counts), 
where it can be seen that they span the same wide range in 
environmental richness. The mean clustering amplitude for 
the RLQs is 213±66 (Model 3), and redoing the RLQ fields 
using Model 1 and 2 for the background counts, the mean 
Bgq's become 214±83 and 224±55 Mpc^ '^'^, respectively. We 
thus find that the mean clustering amplitudes for the RLQ 
and the RQQ samples are practically indistinguishable, i.e. 
on average, there is no difference in the galaxy environments 
on 0.5 Mpc scales for the RQQs and the RLQs. 

The mean and the median Bgq's for the two quasar 



samples are summarized in Table |7|. The median Bgq's are 
seen to be generally lower than the mean values, but this is 
not very significant. The intrinsic scatter seems to be smaller 
for the RLQs, and is more comparable to ABgq than for the 
RQQs. A KS test gives a significance level of 0.72 for the 
null hypothesis that the two samples are drawn from the 
same distribution, i.e. the two samples are consistent with 
being drawn from the same distribution. 

This result is contrary to the traditional picture that the 
environments of RLQs and RQQs differ, i.e. that RQQs are 
found in systematically poorer environments than R LQs. I t 
was reported in the similar study by EUingson et al. (1991c) 



that RLQs inhabited poorer environments than RQ Qs a t 
z ~ 0.4-0.5. However, at z ^ 0.2, Fisher et al. ( |l99e| ) 
and McLure & Dunlop (2000) compare samples of RLQs 



and RQQs matched in redshift and luminosity, and find no 
significant differences, hence our result at 0.5 < z < 0.8 
is consistent with these findings at low red shifts . Also at 
higher redshifts, z ~ 1.1, Hutchings et al. (1995) find ev- 



idence for similar environments around RLQs and RQQs, 
compact groups or clusters of star-forming galaxies. It thus 
appears that both RLQs and RQQs at low to intermediate 
redshifts on average occupy the same type of environment, 
consistent with poorer galaxy clusters, and that the environ- 
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Figure 7. RQQ environment as a function of redshift and Mg . The horizontal line indicat es the amplitude 
of the local galaxy-galaxy correlation function, Bgg = 67.5 Mpc^ '^'^ (Davis & Peebles 1983). 



ments of, at least po-werful, RLQs and RQQs do not change 
with redshift, possibly not even out to z ~ 1. 

Our result is also consistent with the findings that at 
low redshifts powerful RLQs and RQQs appear to be hosted 
by luminous galaxies with luminosities above the break in 
the luminosity function (e.g. Dunlop et al. 1993; McLeod & 
Rieke 1994; Taylor et al. 1996; Bahcall et al. 1997; Boyce et 
al. 1998). Luminous early- type galaxies are known to prefer 
rich environments, and it thus follows that powerful quasars 
hosted by these galaxies also appear in rich environments. 
There are some evidence that the hosts of powerful RQQs 
are massive ellipticals, bu t the results are still somewhat in- 
conclusive. Bahcall et al. (1997) find more RQQs in elliptical 
hosts than in spiral hosts, but note that the quasars also ap- 
pear in_interacting systems and spiral galaxies. McLure et 
al. ( |1999[ ) find that aU Mr < -24 RQQs in their sample 
are hosted by massive ellipticals and Hughes et al. (200C) 



present off- nuclear spectra of luminous AGN, and identify 
the 4000 A break in the stellar component of the host in 
many cases. 

For completeness, we mention that if we exclude the 
three quasars with radio powers > 1.3 x 10^^ WHz~^ (the 



two BQS quasars and LBQS 2348-1-0210), the mean becomes 
191±73 Mpc^'^^ (Model 3), so there is little change, and the 
mean is still consistent with the mean Bgq for the RLQs. 
This is also the case when we exclude the two BQS quasars 
on the basis that the BQS suffers from incompleteness, the 
mean being 177±71 Mpc^'^^ The mean amplitude for the 
sample when the seven radio-detected quasars are excluded 
is 178±82 Mpc^-". 



8 DISCUSSION 

In Paper I, we found that the clustering amplitudes in RLQ 
fields seemed to correlate with radio luminosity instead of 
redshift. Extrapolating this to the RQQs with radio lumi- 
nosities two-three orders of magnitude lower, one might 
expect them to be sited in poorer environments than the 
RLQs. Instead we find that, on average, the environments 
around RLQs and RQQs are indistinguishable. This sug- 
gests that the process that determines the radio loudness is 
not dependent on the environment on Mpc scales, but may 
instead be found in the central regions of the host galaxy. 
This is consistent with the correlation found by Rawlings & 
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Table 6. Mean and median clustering amplitudes for 21 RLQs and 20 RQQs. Model 1-3 correspond to 
different ways of correcting for background galaxies. The intrinsic scatter in the Bgq distribution is denoted 
(Ti„t and was calculated as (cr?, + AB^„)^^^ , where cth^ 



Sample Mean 



Model 1 
Median 



Model 2 Model 3 

Mean Median (Tint Mean 



Median 



CTint 



RQQ 336±77 229 387±54 212±74 151 276±72 210±82 131 315±73 

RLQ 214±83 264 332±75 224±55 165 172±86 213±66 178 236±76 
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Saunders (1991) between narrow-line luminosity and radio 



luminosity in radio galaxies, and the radio-optical correla- 
tion in steep-spectrum RLQs found by Serjeant et al. (1998). 
These correlations link the production of optical emission 
lines in the nuclear regions to the production of large-scale 
radio emission, implying that the fuelling of radio jets is 
linked to accretion onto the black hole. 

If, e.g. as a result of galaxy interactions and mergers, a 
RLQ is born, the environment into which the sour ce expands 
may af fect its radio luminosity, e.g. enhancing it (Barthel & 



Arnauc 199C), and we observe a weak correlation between 



Bgq and radio luminosity (Paper I), but the presence of the 
radio loudness itself would be independent of environment. 
If a RQQ is born, we would perhaps expect a correlation 



between quasar optical luminosity and Bgq (also for RLQs), 
but this would be harder to detect, both because the nuclear 
optical luminosities probe much smaller regions than the 
radio luminosity does, and because the optical luminosity 
may vary. 

We expect a relation between environmental richness 
and quasar optical luminosity since the luminosity created 
by accretion onto the black hole scales with the black hole 
mass. According to the K ormendy & Richstone (1995) and 



the Magorrian et al. (1998 ) relations, the mass of a black hole 



scales with the luminosity (or the mass) of the spheroidal 
component of the host galaxy. (See also Ferrarese & Merritt 
(2000) who find a tight relation between black hole mass 
and velocity dispersions of bulges.) Early-type galaxies, be- 
ing bulge-dominated, and likely to reside in rich environ- 
ments, are therefore expected to host the most massive black 
hole s, and thus the most luminous quasars. McLure & Dun- 
lop ( 2000 ) looked at Bgq vs nuclear luminosity as found from 
the Magorrian et al. relation, but found that the correlation 
was not statistically significant. 

Being sited in equally rich environments, RLQs and 
RQQs are equally replenished by fuel from their surround- 
ings, so both quasar populations would survive for similar 
amounts of time. This implies that the fraction of radio-loud 
quasars would be independent of redshift, and agrees well 
with Goldschmidt et al. ( 199£ ) who find that the radio- loud 
fraction depends on optical luminosity, not redshift. 

There does not appear to be evidence for an epoch de- 
pendence in Bgq between z ~ 0.2 and ~ 0.8. It is however 
not unthinkable that there could be one for RQQs, but at 
higher redshifts, z ~ 1-2. This is because hosts of powerful 
RQQs change with redshift; at low-z, they are giant ellipti- 
cal galaxies with L > L* (the less powerful ones have spiral 
hosts), whereas at high redshifts they are ~ L* galaxies (e.g. 
Hutchings 1995b; Ridgway et al. 2000; Kukula et al. 2000). 
Radio galaxies do not change over the whole interval 2 ~ 2 
to (Lacy, Bunker & Ridgway 2000), so according to radio- 
loud unified schemes, RLQ hosts should behave in th e same 
manner. Aragon-Salamanca, Ellis & O'Brien ( 1996 ) do in 
fact find a significant excess around RLQs at 2 ~ 2.5 that 
is not seen around RQQs at the same redshift (four RLQs 
and six RQQs). 

The change in the hosts of RQQs with redshift is prob- 
ably related to the huge increase in their space density be- 
tween 2 ^ — 1 and 2 ~ 2 — 3. In a model for the form ation 
and evolution of quasars, Kauffmann & Haehnelt (2000) link 



their evolution to the hierarchical build-up of galaxies, and 
predict that quasars at a given luminosity should be found in 
progressively less luminous host galaxies towards higher red- 
shifts. But this explanation cannot work for the radio-louds. 
The fact that the hosts of the radio-louds do not appear to 
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evolve strongly may be because there is a strong link be- 
tween the mass of the black hole and the ability to produce 
powerful jets (Lacy, Ridgway & Trentham 2000; Laor 2000). 
Thus only objects whose black holes (and therefore hosts) 
are already massive at high redshifts can be radio loud. 

In principle, we can also quote a lower limit to the co- 
moving space density of Abell 1-2 clusters at z ~ 0.7. The 
s pace density of q uasars at z = 0.7 is ~ 3 x 10~^ Mpc~^^ 
( Boyle et al. 200C| ) , and since « 10 per cent of the quasars in 
our sample live in Abell 1-2 clusters, this would correspond 
to >3xl0"'' Abell 1-2 clusters per Mpc^ at z ~ 0.7. This is 
potentially interesting for cosmology and may even constrain 
flo, because predictions of cluster nu mber densit ies at high 
redshifts are highly model dependent ( Lilje 199^ . 

It is also interesting to note that this can be used to 
set a lower limit to the duty cycle of quasars. Normally, 
only upper limits are given based on e.g. central black hole 
masses in Seyferts. The duty cycle of quasars can be thought 
of as the probability for a black hole to accrete, or the ratio 
of active to dormant + active quasars at a certain epoch. 
We know from our study that the fraction of active quasars 
living in Abell 1-2 clusters at z ~ 0.7 is « 0.10, and we 
also know that the number of clusters hosting a black hole 
(both dormant and active) must be less than the number 
of Abell 1-2 clusters at the present epoch. If we assume 
that the typical X-ray luminosity of such clu sters is 10"*^ 
ergs ~^, the local cluster density is ~ 2xl0~^ (Rosati et al. 
1998)7jind thus the lower limit to the quasar duty cycle 



is « 0.1x3xl0"''/2xl0"'', i.e. about 1 per cent. This only 
applies to RQQs living in Abell 1-2 clusters, and the present- 
day duty cycle of such systems is probably very low. The 
continous accretion phase cannot last very long, and at some 
point it stops, perhaps when the cluster virializes. 



9 CONCLUSIONS 

The RQQs studied in this paper live in a diversity of envi- 
ronments, from field-like environments to galaxy clusters as 
rich as Abell class 1-2. On average, the sample has a clus- 
tering amplitude corresponding to groups or poorer clusters 
of Abell class « 0, and is indistinguisable from the mean 
amplitude found for a sample of RLQs matched in redshift 
and AGN luminosity. 

Our finding that the average environment around RLQs 
and RQQs at 0.5 < z < 0.8 is statistically indistinguishable 
agrees with wha t is f ound for powerful quasa rs at z ~ 0.2 
by Fisher et al. ([l996|) and McLure & Dunlop (|2000|). There 



is agreement also in terms of the amount of clustering, sug- 
gesting that quasar environments, for both RLQs and RQQs, 
do not undergo much evolution between z ~ 0.2 and ~ 0.8. 
These results are consistent with host galaxy studies finding 
that both powerful RLQs and RQQs live in massive, ellip- 
tical galaxies (e.g. Dunlop et al. 1993; Bahcall et al. 1997; 
McLure et al. 1999) . We therefore believe that the consensus 



i Found by integrating Boyle et al's ( [200o| ) quasar luminosity 
function between the absolute magnitude limits of our survey. 
We used the expression ^pf^^ = [lQ-OA(<^-^){Mu-Ml{z))] ^ 

10-0-4[(/3-i)(A'/b-m* ^hich is different from that given 

in the paper, since we suspect that a misprint may have occurred. 



that RQQs occupy systematically poorer environments than 
RLQs should be put under doubt. 

Since our RQQs on average appear in cluster environ- 
ments they are most likely biased with respect to galaxies. 
Galaxy clusters are also biased (e.g. Lilje & Efstathiou 1988; 
Croft, Dalton & Efstathiou 1999), so the RQQs in clusters 
must have a bias factor that is greater than that for normal 
galaxies. This suggests that RQQs can not be used as unbi- 
ased tracers of galaxies when studying large scale stucture 
using quasars surveys. 

That RLQs and RQQs on average live in equally rich 
environments implies that the origin of the radio loudness 
is not linked to the environments on Mpc scales, but may 
instead be buried within the central parts of the host galaxy. 
This is consistent with the suggestion that the fuelling of 
radio jets is linked to accretion onto a central black hole 
(Rawlins & Saunders 1991; Serjeant et al. 1998). 

A colour-magnitude diagram of the richest RQQ fields 
in this study shows hints of a red sequence at the expected 
colour oi z = 0.7-0.8 early-type galaxies. With spectroscopy, 
these fields will likely turn out to contain galaxy clusters 
a,t z = 0.7-0.8. RLQs and radio galaxies are often used to 
find and study galaxy clusters at higher redshifts, and our 
result that RLQs and RQQs inhabit similar environment 
at z ^ 0.7-0.8 should therefore increase the pool of likely 
cluster candidates at high redshifts. This may be important 
for theories of galaxy formation, since it allows us to find 
and study high redshift galaxy clusters that are not selected 
on the basis of extreme optical or X-ray luminosity. 
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